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We introduce an exactly solvable model of a fermi gas in one dimension and compute the 



momentum distribution exactly. This is based on a generalisation of the ideas of bosoniza- 
tion in one dimension. It is shown that in the RPA limit(the ultra-high density limit) the 
S ' answers we get are the exact answers for a homogeneous fermi gas interacting via a two- 

body repulsive coulomb interaction. Furthermore, the solution may be obtained exactly for 
arbitrary functional forms of the interaction, so long as it is purely repulsive. No lineariza- 
. tion of the bare fermion dispersion is required. We find that for the interaction considered, 

^ ■ the fermi surface is intact for weak repulsion and is destroyed only for sufficiently strong 



0^ 
(N 

o 



repulsion. Comparison with other models like the supersymmetric t-J model with inverse 
square interactions is made. 

Recent years have seen remarkable developments in many-body theory in the form an assortment of 
techniques that may be loosely termed bosonization. The beginnings of these types of tcchiniqucs may be 
\ traced back to the work of Tomonaga 0| and later on by Luttinger and by Lieb and Mattis The 

I I work of Sawada 0| and Arponen and Pajanne in recasting the fermi gas problem in a bose language has 

' to be mentioned. Arponen and Pajanne recover corrections to the Random Phase Approximation (RPA) of 

§ , Bohm and Pines Q in a systematic manner. In the 70's an attempt was made by Luther |Q at generalising 



these ideas to higher dimensions. Closely related to this is work by Sharp et. al. |T^] in current algebra. 
More progress was made by Haldane |^ which culminated in the explicit computation of the single particle 
propagator by Castro-Neto and Fradkin and by Houghton, Marston et.al. |ll| and also by Kopietz et. 
al. Ri gorous work by Frohlich and Marchetti is also along similar lines. Also the work of Frau 

et. al. on algebraic bosonization is relevant to the present Letter as the authors have considered effects 
beyond the linear dispersion in that article. 

In this Letter, we try to use the formalism introduced by us in an earlier preprint to write down 
an interaction jl^], |^ that closely mimics the two-body repulsive interaction between the fcrmions and 
compute the momentum distribution exactly. It is also shown that in the ultra-high density limit, the 
interaction is exactly the two-body repulsive coulomb interaction insofar as it is able to recover the random 
phase approximation exactly. This momentum distribution differs from the traditional Luttinger model in 
one dimension in this aspect, namely, the discontinuity at the fermi momentum is present for weak repulsion 
and is destroyed only for stronger repulsions. But the present model does share similarities with other models 
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such as the supersynimetric t-J model JTzf which does show a discontinuity at the fermi mocmcntum in its 
momentum distribution. 

In the sissa preprint [ p^ we wrote down a correspondence between products of fermi fields and bose fields 
that represent the displacements of the fermi sea, in a manner analogous to the fermi-surface displacements 
introduced by Haldane There, the bose fields in question were written as ak(q) and ak(q) Linear 
combination of these operators would then correspond to displacement oprators of the fermi sea at wavevector 
k by an amount q. This is a generalisation of the concept of fermi-surface displacements where the magnitude 
of k would be restricted to equal kf the fermi momentum. This enables us to draw up a correspondence 
as described in the sissa preprint The boldface on the vectors q and k serve to illustrate that these 

ideas admit a straightforward generalisation to more than one dimension. Let Ck and cj^ be the fermi fields 
in question, then, 



Ck+q/2Ck-q/2 = nF(k)— (5q^o + ( W ) [Afc (q)ak (-q) + Ak(-q)aUq)] 

+ Ak+q/2-qi/2(-qi)Ak-qi/2(q - qi )a^q/2-q^/2 (qi )ak-qi /2 (-q + qi) 

qi 



- Ak-q/2+qi/2(-qi)Ak+qi/2(q - qi )aj,_q/2+qi/2 (^1 )ak+qi /2 (-q + qi) (1) 

qi 

Here, ak(q) and aj^(q) are exact bose operators and Ck and cj^. are exact fermion operators. In symbols, 

{ck,Ck'} = 0; {ck,cj^,} = <5k,k' (2) 

and, 

[ak(q),ak'(q')] = 0; [ak(q), 4' )] = ^K^'^q.q (3) 



and, Ak(q) = ■\/ni?(k + q/2)(l — ni?(k — q/2)) and, ni?(k) — d(kf — k). Furthermore, the kinetic energy 
operator is given in the fermi language as. 



^ = E^k4ck (4) 
k 

where ek = fc^/2m. The same operator has an extremely elegant form in the bose language. 



K = Eo + Y c^k(q)a[(q)ak(q) (5) 

k,q 



here, Eo = ekriF(k), and the excitation energy is given by. 



c^k(q) = (^)Ak(-q) (6) 
m 
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Also the number operator N — cjj.Ck is to be distinguished from its expectation value, (N) = npCk). 
The coefficient Ak(— q) is non-zero only when k.q > 0. The meaning of the above formula in Eq.( |^) is as 
follows. All the dynamical moments of c|j.^q^2'^k+q/2 evaluated in the fermi language are going to be identical 
to those evaluated using the bose representation of this product, provided, we identify the noninteracting 
fermi-sea with the bose vacuum. In other words, 

ak(q)|F5)=0 (7) 

From this it is possible to recast any problem involving interacting fermions in a form involving only bose 
fields. In particular, if we selectively choose parts of the interaction terms written out in the bose language, 
it may be possible to compute exactly, important quantities such as, momentum distribution, dynamical 
density correlation functions and so on. Let us now focus on the homogeneous fermi gas. The interaction 
has the form, 

q/O 

Written out in terms of the bose fields, it has a part that is quadratic in the bose fields, and we have cubic 
and quartic terms as well. Let us just focus on the part that is quadratic in the bose fields. This amounts 
to postulating a phenomenological interaction that is not exactly of the the two-body type (in the original 
fermi system) but something that mimics it. How closely does it mimic this two-body interaction, is the big 
question. This question may be answered by computing the dielectric function and demonstrating that it is 
exactly equal to the RPA dielctric function. This has been done in our important preprint that computes 
the single-particle green functions of the same system Thus our results are the exact results for a 

homogeneous fermi gas with two-body repulsive interactions in the same limit in which the RPA is exact. 
Written out in the bose language this interaction has the form, 

= E 1^ 5][Ak(q)ak(-q) + Ak(-q)a[(q)] [A^' {-ci)a^' (q) + A^, {ci)al, (-q)] (9) 
q#o k,k' 

The full hamiltonian is now H = K + Hj . This hamiltonian may be diagonalised exactly and the momentum 
distribution may be computed. Alternatively, the equation of motion method may be used to obtain the 
same answers. In fact, we have used both these methods po| to obtain the exact answer for the momentum 
distribution in one dimension (the formulas written down till now are valid in any number of dimensions). 
The exact momentum distribution found by us has the form {h = c = I), 

dqi Ak-q^/2{-ql) 



211 2ujR{q,){ujR{q,) + Uu-q,,2{qi)?{^){cosh{\{qi)) - 1) 



{2'Kk ) ^ ^k+q^/2{-qi) 

' , 27: 2un{qi){LUniqi)+u;k+q,/2{qi))H'^)icosh{X{qi))^l) 



m = {—)(-) (11) 

m Vq 



Jgk (fc/ + - (fc/ - g/2)2exp(-A(g)) 
-i^.('Z) = (-)y l-exp(-A(g)) (^2) 

The above answer for the momentum distribution Eq. ( ^) is the exact momentum distribution of a homo- 
geneous fermi gas in one dimension interacting via a two-body repulsive interaction. As has been pointed 
out earher, the RPA dielctric function as descibed by Bohm and Pines may be recovered exactly by this 
model provided one identifies the density operator with the object obtained by retaing only terms linear in 
the sea-displacements. 

Pq = ^[Ak(q)ak(-q) + Ak(-q)a|,(q)] (13) 
k 

The above form ignores the terms quadratic in the sea-displacements. The justification being that this 
definition along the with the standard procedure for obtaining the dielectric function, namely, add a periodic 
weak potential coupled with the density operator above and compute the ratio of the external versus the 
effective potential, gives the exact RPA dielctric function thereby suggesting that dropping terms quadratic 
in the displacements (in the definition of the density) is not going to be important in the same limit in which 
the RPA itself is exact (for a pedagogial description of this procedure of obtaining the dielectric function 
please consult the text by Kadanoff and Baym [Q). In order to see how good the present theory is, it is 
desirable to compare these results with the Calogero-Sutherland model or more specifically with the spin- 
spin correlation function of the Haldane-Shastry model. This is given by Lesage et. al. Here they 
find that the momentum distribution is smooth and does not possess any discontinuity or kinks. This is 
in contrast with the solution above (Eq.( ^^), where we see the presence of a fermi surface for sufficiently 
weak repuslsion. However, the supersymmetric t-J model |T^ ] with inverse square interactions does posses a 
jump at the fermi surface, and therefore the present model belongs to this latter class. Also, the system does 
not exhibit Wigner crystallization at low density. This is another drawback of the interaction that we have 
considered. But this is to be expected since the present model captures the two-body repulsive interaction 
only in the ultra-high density limit. On the plus side, the highly nonperturbative nature of the solution 
suggested by the presence of the term, cosh{^^-^) — 1 cannot be missed. 
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